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Abstract 

From a simple path integral involving a variable volatility in the veloc- 
ity differences, we obtain velocity probability density functions with exponen- 
tial tails, resembling those observed in fully developed turbulence. The model 
yields realistic scaling exponents and structure functions satisfying extended 
self-similarity. But there is an additional small scale dependence for quantities 
in the inertial range, which is linked to a slow approach to Kolmogorov (1941) 
scaling occurring in the large distance limit. 



The universal features displayed by fully developed hydrodynamic turbulence are 
still not fully understood. Kolmogorov (1941) |I|, showed how a set of statistical 
quantities known as structure functions are expected depend on the length scale r 
as power laws with predicted exponents. Experimental measurements have indicated 
that while these predictions are close to the truth, the predicted exponents are not 
exactly realized. In the face of this experimental input, much effort has been devoted 
towards understanding the origin of these anomalies in the scaling exponents, while 
retaining the notion that current experiments are observing an "inertial range" where 
strict power law scaling holds. 

In this paper we will investigate what appears to be a loop-hole in this reasoning. 
In spite of impressive advances made in the experimental studies, the fact remains 
that the scaling exponents have been deduced by looking at scaling regions extending 
over little more than one decade in r || [|, ||, [(J. This leaves open the possibility 
that a small but significant departure from strict power law scaling is still consistent 
with the data. We will argue that this possible departure is sufficient to allow the 
observed anomalies in the exponents to be nothing more than a transient effect related 
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to intermittency, and that the true power law scaling only occurs on larger scales. 
This very large distance scaling could take the form proposed by Kolmogorov. 

We shall present a model which predicts deviations from power law scaling, and 
which shows that these deviations may be small enough to have escaped detection 
thus far. The model is based on a simple physical picture for the effects of inter- 
mittency, which are effects due to the coming and going of coherent structures in 
the velocity field. We will model the effects that these structures (eddies) have di- 
rectly on the probability distribution function (PDF) for velocity differences. We are 
claiming simplicity in our approach, but not uniqueness. Thus the applicability to 3- 
dimensional hydrodynamical turbulence in particular could be considered speculative, 
since the model makes no mathematical contact with the Navier-Stokes equations. 
On the other hand, when the consequences of the model are explicitly worked out 
they are found to coincide rather well with the data. 

Only the symmetric part of the observed PDF will be modeled in this paper. This 
is denoted by P r (Sv r ) where 5v r is the difference in some velocity component at two 
points separated by a distance r. The model will yield an explicit expression for the 
PDF which 1) for small r has the typical sharp peak and broad tails characteristic of 
intermittent behavior, 2) at any finite r displays exponential tails at large enough 5v r , 
and 3) tends toward a Gaussian form in the large r limit. But the model shows that 
these realistic features may in fact be implying that the values of the scaling exponents 
are slightly scale dependent, as they evolve from their observed "anomalous" values 
for the values of r where they are presently measured, to the pure K41 values in the 
large r limit. 

We consider a discrete set of points on the line connecting the two points separated 
by distance r. We label the points by their distance from one end, (0, n, r 2 , rjv, rw+i) 
where rjy+i = r. We start by assuming that P r (5v r ) can be approximated by 
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with Pf^ V(y)dy = 1. We assume that the set of points rj is characterized by a scale 
p such that {rj/p) a = j, where j = 1, 2, N + 1 and a is a positive constant to be 
determined below. The integration over all velocities at the intermediate points may 
be thought of as a sum over all paths in the space of velocities which start at vq and 
end at v r . We could thus refer to ([!]) as a path integral. But note that we refrain from 
taking the continuum limit, N —>■ oo with r fixed, since the scale p has a physical 
significance. 

The discretization of the range from to r into subregions (rj, Tj + i) is associated 
with our modeling of intermittent behavior. Roughly speaking, we are suggesting 
that the coming and going of eddies, of sizes comparable to subregion size r J+ i — rj, 
produce a variable volatility in the typical velocity differences v rj+1 —v rj . We make this 
statement precise by writing V(v r . — i> r -_i) as a superposition of Gaussians, where we 
integrate over all values of a volatility parameter Sj which is itself weighted according 
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to a Gaussian. To simplify notation we denote v r . 
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The consideration of large and small values of Sj corresponds to the possible presence 
or absence of eddies in that subregion at various times. Note that the variance of 
both V(xj — Xj-i) and -PGauss(0", x) is a 2 . 

([I]) is constructed to give the cumulative effect of the variable volatilities occurring 
in the various subregions. p should be typical of scales on which viscosity influences 
and damps the formation of eddies. As such p is expected to lie between the dissipative 
(Kolmogorov) scale, rj, and the lower end of the inertial scaling range (the latter range 
being characterized by the absence of viscosity effects). Note that if < a < 1 then 
the size of the subregion, r J+1 — rj, and thus the relevant eddy size for that subregion, 
would increase with j. 

We proceed by inserting (0) into (0) and integrating over the Xj. We can write 
the terms in the exponential which depend on y = (x\,X2, ...,x^) in matrix notation 
y T My + Jy where M is a matrix and J is a vector, with the latter depending on xq 
and xn+i- We may complete the square and do the Gaussian integrations over y, 
for various values of N. From this we are able to deduce a simple result, where all 
dependence on the Sj now resides in S 2 = Ylf=i s "j- The Sj integrations for fixed S 
then just give the surface area of an iV + 1 dimensional sphere. 



P N (x N+ i - Xq) 



<t(27t)t+i 

1 

a(27r)T+ 1 \ r 



N+l 

n 

3=1 




(XN+I — XqY 



2cx 2 S 2 




x / dSS"- L exp 



{XN+I — Xq) 2 _ 1^2 



2a 2 S 2 



cr0Fr(^±i) h 



{x N+ i - x y _ i 

2a 2 ^ 2 (iV + l) 2 



S 2 (N + 1) 



(4) 



In the last step we have defined S 2 = S 2 / (N + 1) for convenience. 

Finally we will take this result and analytically continue from integer values N + l 
to positive real values r. Replacing iV + 1 by r = (r/p) a givesQ 
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lr This superposition of Gaussians is reminiscent of the proposal in Q|. 
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We note that for large r the integral over S becomes strongly peaked about S = 1, 
and thus 

P r (5v r ) -> — 7 =exp • 6 

In addition we find, for any r, that the variance in the velocity differences is 

({5v r y)=o 2 T = o-\r/pY. (7) 

The variable r controls the evolution of the PDF, in the sense that evolution 
through unit steps in r are generated by successive application of the "evolution 
operator" in (0). In this sense r is an "evolutionary" time scale depending on the 
distance scale r. There is also a dynamical time scale associated with scale r, the 
eddy turnover time, which is simply r/v where v is a typical velocity observed on 
scale r. If we associate v with G\Jt from (|7|), and if the evolutionary and dynamical 
times happen to scale with r in the same way, then we obtain a = 2/3. This gives 
the standard scaling law ((Sv r ) 2 ) oc r 2//3 . 

We now perform the integral over S, and rescale the PDF such that the variance 
is unity for all r. We obtain 

PAX) = > l' ( ^ 1) *i(T-l)(HT l/a ). (8) 
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K v (y) is the modified Bessel function of the second kind. In Fig. (1) we display P T (x) 
for various r. 

When t is an even integer, P T (x) is an exponential times a polynomial in x, and 
in particular for r = 2 it is purely exponential, /^(x) = exp(— \x \V2)/y/2. For 
t < 2 the PDF is even more strongly peaked and it decreases less quickly than an 
exponential ("stretched exponential"), before becoming purely exponential at large 
enough x. For r > 2 we have "deformed Gaussians" with exponential tails. For 
any r, dlnP T (x)/dx —y/r. Alternatively we may consider moments, M™ 

j P r {x) \x\ n dx, which are all finite and given by 

on -p ( n+1 \-p ( n+r \ 

m; = * 1 2 ; 2 } . (9) 

The approach to the exponential tails is reflected in M^ +1 /((n + 1)M") n z^° 1/^/7. 

It is common to describe the evolution of the PDF with r in terms of the gener- 
alized structure functions, defined by S n (r) = (\5v r \ n ). For r in the inertial range it 
is usually assumed that S n (r) oc r^ n . Our model suggests that exactly scale indepen- 
dent exponents do not exist for the typical ranges of r considered, and that instead 
we should consider the local exponents ( n ( r ) = d\og(S n (r)) / ^log(r) . We obtain 

<n(r) = |r(«r(i(n + r))-*(ir)), (10) 
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Figure 1: log 10 [P T (x)] vs. x for r = 1/2,1,2,4,8,16,32, top to bottom. The curves 
are vertically displaced by 0.5, 0, —0.5, —1, —1.5, —2, —2.5 respectively. 



where \1/ is the digamma function. This gives C2( r ) = a i and in the large r limit, 
£ n (oo) = an/2. The latter is K41 scaling when a = 2/3. Although we adopt this 
value of a in the following, we leave open the question of whether a is exactly equal to 
2/3 or just close to it.0 We also stress that we are describing only the symmetric part 
of the observed PDF, which in the case of the PDF for longitudinal velocity differences 
also has an asymmetric part. This asymmetry is reflected by the nonvanishing of the 
structure functions (Sv^ for odd n. In particular our £ 3 is not the same as the 

scaling exponent of (8v^\ which is constrained to be unity 0, |^]. 

In Fig. (2) we display Cn(^) for various r, and compare to the large r limit values 
Cn = n/3. The main point is that the scaling exponents approach their asymptotic 
values very slowly. It is also of interest that for some range of r/p, our exponents are 
realistic for a wide range of n. We show this by comparing to the model of She and 
Leveque (SL)@, which yields = n/9 + 2 — 2 (2/3) n ^ 3 , and which is known to fit the 
current data quite well. We may phrase the agreement in terms of the relative scaling 
exponents Cn/C3> which are known to better accuracy [||, f|, ]5|, || than the individual 
exponents, and for which the factor a in ( |I~0"| ) cancels. We find, for example, that 
Cn(35p)/Cs(35p) as a function of n is within a few percent of the corresponding SL 
values up to n = 30! And as r varies from 20 to 50, Ce( r )/C3( r ) f° r example varies 
from 1.72 to 1.82, compared to C 6 Si /C 3 SL = 1-78. 

To see the r dependence of the structure functions themselves we display S n (r) in 



2 For example the first reference in |3| advocates (2 = .7. 
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Figure 2: ( n (r) vs. n for r/p = 10,20,40,80, 160, bottom to top. The dashed line is 
Cn = n/3. 



Fig. (3). Over some range of r the lines appear to be close to straight. But other than 
for n = 2, which is exactly straight, the lines have a positive curvature which increases 
for larger n. We will see below how these positive curvatures, or equivalently the scale 
dependence of the scaling exponents, may be rather subtle to observe experimentally. 
But first we must consider the effects of viscosity on the structure functions. 

We expect that the effects of viscosity will cause r = (r/p) 2 / 3 to be replaced by 
a function /(r) which deviates from (r/p) 2 / 3 as r approaches the dissipative scale 77 
from above. Since a decreasing r corresponds to increasing intermittency, and since 
the cause of intermittency in our picture — variable volatility of velocity differences — 
should be damped by viscous effects, we expect that /(r) should start to decrease 
more slowly than (r/p) 2 / 3 . Thus the evolution of the PDF shapes is retarded by 
viscous effects in the "intermediate viscous range" , the range of r between 77 and the 
start of the inertial range. The implication is that at large scales in the inertial range 
where f(r) = (r/p) 2//3 it is likely that p > 77. We may also expect that p/77 grows 
with the size of the intermediate viscous range. 

For very small r where the velocity field is smooth the variance is evolving like r 2 , 
faster than in the inertial range, and thus there must clearly be a decoupling between 
the evolution of the PDF shapes and the evolution of the PDF variance. We are thus 
led to a generalization of our model, where we replace the expression in ([|) by the 
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Figure 3: log 10 [S' n (r)] vs. log 10 (r) with p = 1 for n = 2,4,6,8,10, bottom to top, 
each line vertically offset so as to vanish at r = 1. 



following.^ 

p <Sv ) - (f(r)m' J(T, - m r ass'V-'exo ( - Is° fir)) ill) 

This new PDF has the property that ((5v r ) 2 ) = g(r), and thus we can use g(r) to 
reproduce the observed behavior of the variance even in regions where viscosity or 
finite size effects are important. The function f(r) determines the evolution of the 
PDF shapes; that is, the new P r (x) is obtained from the one in (|5p by replacing r 
with /(r). 

We will illustrate the effects of viscosity with specific choices for f(r) and g(r). We 
consider /(r) = ([r + p — g]/p) 2 ^ 3 , which is unity at r = r\ and approaches (r/p) 2 / 3 
at large r.Q For g(r) in the range rj < r < 10 A rj we consider 



9(r) 



(r/77) 



2/3 



3 This expression could be "derived" from a path integral as before, where the discretization of 
the range from to r would be determined by unit steps in the function g(r). Instead of each point, 
groups of adjacent points would have a single fluctuating volatility, and the transition from one 
grouping to the next would determined by unit steps in the function f(r). 

4 Although we expect f(rf) to have some dependence on the size of the inertial range, we take 
/ (77) = 1 for simplicity. 
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Figure 4: \og 10 [S n (r) / r^ 35 ^] vs. log 10 (r) with 77 = 1 for n = 3,6,8,10, bottom to 
top. The solid and dashed lines correspond to g\ (r) and g2^r) respectively. 



We take two examples for (01,02,03,0^,05) which idealize typical data sets, gi(r) 
with (80, 79, 60, 10, 5) having a relatively large inertial range (large Reynolds number) 
and g2{r) with (80,79,70,5,30) having a small inertial range. Our results are not 
very sensitive to precisely how these functions deviate from r 2 / 3 behavior outside the 
inertial range. Note that the intermediate viscous range is larger for gi(r)|] Due to 
our expectation that p should increase with the size of the intermediate viscous range 
it is natural that p\ > P2. For illustrative purposes only, we choose p\ = 4.3r/ and 
P2 = l-3?7. 



With this input we can extract all higher order structure functions from (|TT|) . The 



scaling in the inertial range for both cases turns out (because of our fortuitous choice 
of pi and P2) to be described by the realistic scaling exponents ( n (35p) from (|I0|) . We 
show this in Fig. (4) where we plot S n (r)/r^ n ^ 35p ^ .[] Of most interest are the n = 8, 10 
curves for the g\(r) case, where we see (very slight!) evidence of a double hump 
structure growing more prominent with increasing n. This positive curvature in the 
middle of the scaling region is arising from the positive curvatures in Fig. (3), which 
in turn is a reflection of the gradual approach to K41 scaling. This is the generic 
signature of our model, which should be seen when the size of the inertial range and 
the order n are both large enough. The figure makes it clear, though, that it is a 

5 The tendency for the intermediate viscous range to grow with the Reynolds number has been 
noted for example in the first reference in ||. 

6 This figure may be compared with Fig. (1) in the second reference of Q. 
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Figure 5: log 10 [5 n (r)] vs. log 10 [S , 3 (r)] for 4 < r/r] < 10 4 and for n = 6,8, 10, bottom to 
top. The lines for the g±(r) and g2^r) cases are displaced by 0.5 and —0.5 respectively. 
The dashed straight lines have slopes ( n (35p)/( 3 (35p). 



small effect which may be hidden in present data. 

Our model may also be used to illustrate "extended self-similarity" JTI], [|. In Fig. 
(5) we plot S n (r) versus S 3 (r) for 4 < r/r] < 10 4 for the two cases g\{r) and g2{r). 
For comparison we add straight lines with slopes given by the relative exponents 
Cn(35p)/Cs(35p). We see that scaling has been extended to smaller scales than is 
apparent in Fig. (4). Similar results are obtained for other choices of the functions 
f(r) and g(r). It thus appears that such plots are not very sensitive to the deviations 
from power law scaling we are proposing. 

We reiterate that the model yields a universal set of local scaling exponents, and 
the scaling exponent from an experiment depends on what distance scale, effectively, 
the local scaling exponent is being measured at relative to p. On the other hand if it 
is true that p/r) increases with the size of the intermediate viscous region, as we are 
suggesting, then it may be difficult to obtain measurements at distances, in units of 
p, which are very different from each other. The variability in the scaling exponents 
is also obscured if different experiments have intermediate viscous regions of similar 
size, and/or if one is confined to the lower order structure functions, such as n = 6 
and below. 

We have seen how viscosity and finite size effects can have the effect of transform- 
ing the structure functions in Fig. (3) into the structure functions in Fig. (4) which 
in turn display extended self-similarity in Fig. (5). It is also encouraging to find that 
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realistic values of the scaling exponents emerge when the p parameter is within the 
intermediate viscous range. But perhaps most important is that the model suggests 
how evidence for a slow approach to K41 scaling could eventually be uncovered. We 
should also differentiate between the evolution of the PDF shapes as a function of 
r (to which the structure functions are sensitive), and the basic set of shapes pre- 
dicted by the model, given in (||) and depicted in Fig. (1). These probability density 
functions may be of interest in various other contexts. 
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